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ABSTRACT 



Aims. We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used 
for a wide range of applications of spectral analysis and stellar parameter synthesis. 

Methods. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra 
from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. 

Results. The synthetic spectra cover the wavelength range from 500 A to 5.5 yum with resolutions of R = 500000 in the optical and 
near IR, R = 100 000 in the IR and A/i = 0. 1 A in the UV. The parameter space covers 2 300 K < r cff < 12 000 K, 0.0 < log g < +6.0, 
-4.0 < [Fe/H] < +1.0, and -0.2 < [a/Fe] < +1.2. The library is a work in progress and we expect to extend it up to T c ff = 25 000 K. 

Key words. Convection - Radiative transfer - Stars: abundances - Stars: atmospheres 



1. Introduction 

Stellar spectral libraries are an important tool for analysing 
observed spectra and can be used for a wide range of ap- 
plications. The use of empirical libraries guarantees real 
physical spectra, which will match other observations well, 
e.g. absorption lines created by atomic and molecular transi- 
tions appear at the exact position with the correct shapes and 
strengths. However, these spectra are limited both in spectral 
resolution and wavelength coverage - and extensions of the 
library require some significant effort and observing time, if 
possible at all (e.g. changes in instrumentation). Furthermore, 
the chemical composition of the observed stars is not known 
exactly, resulting in possible systematic errors for abundance 
determinations, and the spectra cover only a limited parameter 
space, usually restricted to the effective temperature, the surface 
gravity, and the metallicity. Examples for empirical libraries 
covering the optical wavelength ra nge are the Indo-U. S. Library 
of Coude Feed Stellar S pectra dValdes et all l2004 MILES 



(fl97l . from 
rived (e.g. 



which many spectral librari es 
Rodriguez-Merino et al.l (|2005l) it 



have 
i the 



or uouae feed stellar spectra (Vaiaesetai. zuiM j), miles 
dSanchez-BTazquez et all |2006|). STELIB (iLe Borgne et alj 
[2003b, and ELODIE dPrugniel & Sojjbjr^|2004. Note thatliil 
these libraries contain only low to medium resolution spectra. 

Synthetic spectra on the other hand are limited by the 
completeness of the spectral line lists and by the knowl- 
edge of line broadening parameters and numerical assump- 
tions - e.g. plane-parallel versus spherical geometry and lo- 
cal thermodynamic equilibrium (LTE) versus statistical equi- 
librium (NLTE). Nevertheless, synthetic libraries have many 
advantages: the range of stellar parameters, elemental abun- 
dances, and both the wavelength range and the spectral res- 
olution can be adjusted as needed. The best know n set of 
model atmospheres currently available was created by iKuruczl 



been de- 
UV, and 

iMurphv & Meiksinl (I2004 and iMunari et al.l (12005I) in the full 
optical domain ). The NLTE stellar atmosphere code TLUSTY 
(lHubenvll988[) made possib le a grid of NLTE a t mosph eres for O 
stars (ILanz & Hubenvll2003l) . lGustafsson et al.l(l2008l) published 
a grid of spectra based on MARCS which, like our new grid, 
uses spherical-symmetry for its models. 

Because we need a stellar s pectral library for t he analysis 
of spectra observed with MUSE (Bacon et al. 2010), X-shooter 
dVernetet al.ll201ll) . and CRIRES dKaeufl et alJl2004h for exam- 
ple, a decision was made in favour of the flexibility of a syn- 
thetic library that can be extended on the fly - even in param- 
eters that are not changed in the published grid. For example, 
varying the abundances of single elements enables us to deter- 
mine the areas of the spectra that are mainly affected by those 
elements. In comparison with empirical spectra, this shows a 
major advantage of having a stellar atmosphere code available 
together with thousands of synthetic atmospheres that can be 
used as starting points for extensions. Instead of extending an 
existing synthetic library, we opted for creating a new one using 
the latest st ate-of-the-art version of the stellar atmosphere code 
PHOENIX dHauschildt & Barodll999h . 

The first published PHOENIX libra ry consisted of 700 
mode l atmospheres for M (sub)dwarf stars (Alla rd & HauschiTdl 
1 19951) . all calculated under the assumption of local thermody- 
namic e quilibrium. This has also been the case for the NextGen 
library (Hausc hildt et al.iri999al) with its temperatures ranging 
from 3000 K to 10000K, which has been used as the starting 
point for a grid of NLTE model atmospheres of dwarfs and gi- 
ants. Furthermore, an upgrade to spherically symmetric models 
for giant stars with effective temperatures between 3000 K and 
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6800 K has been done bv lHauschildt et alJ dl999bh. Further cal- 
culations of NLTE models have been conducted bv lBarman et al.l 
(2000). The effects of dust settli ng in cool s t ars lik e brown 
dwarfs have been investigated b y lAllard et al.l (2001) and ex- 
tended to extra-solar planets by Witte et a ll (120091 1201 lj). The 
latest comprehensive grid of spectra dBrott & Hauschildt 2005) 
was c reated for use with the GAIA mission (Lindegren et al.l 
120081 In the following we will refer to spectra created with that 
version of PHOENIX as AMES-cond-v2.6. 

Our new spectral library is the first one using the current 
version 16 of PHOENIX, which has a new equation of state as 
well as an up-to-date atomic and molecular line list. This al- 
lowed us to produce spectra that match observations especially 
of cool stars significantly better than other synthetic libraries. A 
full dis cussion of this aspect will appear i n a subsequent publi- 
cation dWende - von Berg & Reinersl2013l) . The use of spherical 
geometry guarantees a consistent model grid from the main se- 
quence up to giants. We have adjusted the wavelength range and 
resolution to meet the requirements for science cases using ex- 
isting and upcoming instruments. 

In this paper we will first describe the basic parameters of 
the grid, namely its dimensions and the resolution of the spec- 
tra. Then we will discuss in detail some of the input parameters 
such as the used element abundances, the stellar mass, and the 
micro-turbulence for which we will also provide comparisons 
with observed data. We will characterise the equation of state 
used for calculating the models and show the improvements that 
we achieved compared to older PHOENIXmodels, other spec- 
tral libraries, and observations. Finally, we will explain how to 
download the FITS files containing the spectra and describe their 
content. 

2. The library 

2. 1. Dimensions and resolution 

Since the library is still a work in progress, we only present its 
current state which, however, already covers the parameter space 
of most stellar populations that are not part of starburst regions. 
Table Q] shows the current parameter space of the grid, a fu- 
ture extension will probably be towards higher effective temper- 
atures, most probably up to 25 000 K. The grid is complete in its 
first three dimensions T e e, logg [cgs] and [Fe/H]. So far, a vari- 
ation of the alpha element abundance {[a/Fe] + 0) has only been 
calculated for effective temperatures 3 500K < < 8 000K 
and metallicities -3 < [Fe/H] < 0. 

The standard spectroscopic abundance notation is used in 
this article which defines the ratio of two elements A and B in 
a star relative to their ratio in the Sun as 

[A IB] = log[»(A)/»(B)] - log[n(A)/n(B)] (3 , (1) 

with the number density n. We assume that [X/H] = [Fe/H] for 
most elements and so, [Fe/H] denotes the overall metallicity. 
The abundances of the a elements (O, Ne, Mg, Si, S, Ar, Ca, 
and Ti) are defined as 

[a/H] = [Fe/H] + [a/Fe]. (2) 

This implies that the value of [Fe/H] is conserved in alpha- 
enhanced or depleted models, but the overall metallicity Z 
changes. 

The spectra cover the whole wavelength range from 500 A 
up to 5.5 yum; the sampling is given in Table [2] Please note that 
these sampling rates are the direct output of PHOENIX, i.e. the 
spectra have never been resampled or convolved with any kernel. 



Table 1. Parameter space of the grid. Alpha element abundances 
[a/Fe] + are only available for 3 500 K < T e $ < 8 000 K and 
-3 < [Fe/H] < 0. 



Variable 


Range 


Step size 


7* [K] 


2 300 - 7 000 


100 




7 000- 12000 


200 


logg 


0.0 - +6.0 


0.5 


[Fe/H] 


-4.0 - -2.0 


1.0 




-2.0- +1.0 


0.5 


[a/Fe] 


-0.2- +1.2 


0.2 



Table 2. Sampling of the spectra in the grid. 



Range [A] Sampling 
500 -3 000 A^ = 0.lA 
3 000 -25 000 500000 
25 000 -55 000 R~ 100000 



The wavelength grid is identical for all models in the grid, 
which was made possible by a new option in PHOENIX vl6, 
which has been introduced in this new version and allows us to 
disable the creation of new wavelength points to better account 
for opacity. We activated this option, since the sampling rate is 
high enough already. Furthermore, it allows us to store the wave- 
length grid only once for all spectra in the library and therefore 
save storage space. 

2.2. Calculation of model atmospheres and synthetic spectra 

Synthesising a spectrum with PHOENIX is a three-step process. 
First, a model atmosphere was calculated using a set of 766 
atmospheres from a previously unpublished PHOENIX grid as 
start values. From those, the grid was extended step by step, al- 
ways using an existing neighbouring model as the starting point. 
The PHOENIX code iterates a model atmosphere, until its cri- 
teria for convergence are reached. In order to limit the run time, 
we stopped after every 80 iterations and restarted PHOENIX if 
necessary. Parameter adjustments for many of the models had to 
be done individually in order to obtain fully converged model 
atmospheres. 

The next step was to activate a mode in PHOENIX that trig- 
gers the use of special line profiles, which we used for the Ca 
lines, for example. In this mode no convergence is checked, so 
we carried out five more iterations by default in order to strat- 
ify the atmosphere again properly. Finally, the high resolution 
spectrum was synthesised from this intermediate model. 

Having access to the model atmospheres as well as the spec- 
tra allows us to quickly expand the grid in any given direction. 
Furthermore, we can use the existing model atmospheres for 
other investigations - e.g. we activated the output of the full ra- 
diative field, which can be used for calculating limb darkening 
coefficients in PHOENIX for some models. The results will be 
discussed in a subsequent publication. 

Even in its current state (end of 2012) with ~50000 model 
atmospheres and corresponding synthetic spectra (of which 
-30 000 will be published), it took about -135 CPU years for the 
new PHOENIX library to be calculated on the Nehalem Cluster 
of the GWDG0 in about one year, parallelised on 20-50 nodes 
with 8 CPUs each. 



1 Gesellschaft fur wissenschaftliche Datenverarbeitung mbH 
Gottingen, http: //www. gwdg . de/index.php?id=2156 
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2.3. PHOENIX settings 



Table 3. Coefficients for calculating the mass in Eq.[6] 



While most other synthetic libraries set many parameters to fixed 
values, we tried to use parametrisations, for the mass and the 
mixing length parameter a, for example, that better match ob- 
servations. Below we introduce a new unique method for the 
micro-turbulence, which we calculate from convection veloci- 
ties obtained from the atmospheric model. Furthermore our new 
stellar library is the first one to use the new equation of state 
introduced with the new version of PHOENIX and we use the 
most recent values for the solar element abundances. 

All model atmospheres have been calculated in a spherical 
mode with 64 layers. The reference wavelength defining the 
mean optical depth grid, which should be set to a wavelength 
with sufficient flux, is fixed to A T = 12 000 A for T eS < 5 000 K 
and A T = 5 000 A for hotter stars. 

Local thermodynamic equilibrium (LTE) has been assumed 
for all the models in the library, which is sufficient for the tem- 
perature range that we have investigated so far - the extension 
to effective temperatures T e ff > 12 000K will take into account 
NLTE effects. Nevertheless, even for the existing models with 
T eff > 4000K, we used the NLTE mode of PHOENIX to use 
special line profiles for some species (Li I, Na I, K I, Ca I, Ca 
II). The full sample of molecules is only used for models with 
T e ff < 8 000 K, while it is stripped to a few important ones for 
hotter stars. 

Condensation is included i n the equation of s tate but ignored 
in the opacity (Cond setup, see lAllard et al.l l2001) and we do not 
handle dust settling at all in our model atmospheres, since all our 
models have effective temperatures r e ff > 2 300 K. 

2.3.1. Mass 

While plane-parallel geometry is a valid assumption for main 
sequence star model atmospheres, this is no longer true for gi- 
ants. For those, PHOENIX offers a mode for spherical symme- 
try, which we used for all models to guarantee a homogeneous 
grid. A model atmosphere with spherical geometry is then de- 
scribed by its effective temperature r e ff , surface gravity log g, 
and mass M*, from which all the other physical parameters can 
be derived, i.e. adjusting the mass changes both the radius and 
the luminosity. In order to get numerically stable and realistic 
atmospheres, we used a varying mass throughout t he gri d and 
did not fix it to specific values (as lHauschildt et al.l d!999bl) and 
Gustafsson et al. (2008) do, for example). 

For main sequence stars we use the mass-lumi nosity rela tion, 
which can be approximated for solar type stars (Voigt] [l99Tl) as 



LJL = (MJMof 



(3) 



Furthermore we can calculate the luminosity from the effective 
temperature and radius. For the Sun the radius can be substituted 
by the surface gravity and mass, which yields: 



LJL B = (T eS ,JT eS , B )\RJR Q ) 2 = (T eS ,JT e ^ ) 4 MJM Q . 
By combining these two equations we get: 

(r eft >/r eff , ) 4 = (mjm q ) 2 . 



For evolved stars, we can naturally not prescribe a unique 
mass-radius relation, and it becomes more challenging to define 
mass as a function of T e ff and log g, so we modify the main se- 
quence relation with a variable factor c for lower surface gravi- 
ties to match observed values of giants and super giants. Masses, 
especially of red giants and asymptotic giant branch stars, show 



log(g) 


> 4 > 3 


> 2 


> 1.6 


> 0.9 


> 
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1 1.2 
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Fig. 1. Distribution of stellar masses for different effective tem- 
peratures T e ff and surface gravities log g, colour-coded in units 
of solar mass from M Q (black) to 25 M Q (white). 



considerable spread even for identical atmospheric parameters, 
but we chose adjustment factor s such as matching the very-well- 
constrained values of Arcturus (iRarrurez & Allende Prietol201 ll) 
and the sample of M giants studied by Dumm & S child! ([1998), 
roughly covering 3000 K < r eff < 4000 K, -0.5 < logg < 1.5 
with masses between 0.8 and 4M . Therefore in the grid the 
mass is given by 



M* = c ■ Mq • (r e ff/5 770 K) 2 



(6) 



with values for the coefficient given in Table [3] Figure [T] shows 
the distribution of masses over the grid. For solar-like stars, we 
obtain = 1M and the masses increase for hotter and larger 
stars as expected. 

2.3.2. Convection 

In PHOENIX, convection in the stellar atmosphere is described 
by the mixing length th eory of turbulent transport, whi ch goes 
back t o iPrandtll d 1925b and, in the stellar context, to IVitensd 
dl953l) and B6hm-V itenseldl958l) . It is characterised by the mix- 
ing length parameter a, which describes the ratio between the 
characteristic length that a volume of gas can rise in a stellar at- 
mosphere before mixing with its surrounding and the pressure 



(4) scale height (a = l/H p ), i.e. it is a value for the efficiency of 
the convective energy transport - a large a indicates an efficient 
energy transport, while it is inefficient for a small a in which 
case the gas volume is allowed to rise to greater heights under 
the buoyant force, thus reaching higher velocities and carrying 
more heat before it is assumed to be dissipated. 

For our sp ectra we descr i be the mixing length parameter a 
as discussed in lLudwig et al.l d 19991) where it has been calibrated 
using 3D radiative hydrodynamic models. Since the authors em- 
phasise that it might be incorrect to extrapolate the mixing length 
outside the ranges given in the paper (4 300 K < T e g < 7 100 K 
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3000 4000 5000 6000 7000 8000 
T eff [K] 



Fig. 2. Distribution of the mixing length parameter a in the grid. 



and 2.54 < logg < 4.74), we restrict the result to an interval 
a e[l,3.5]. 

The distribution of mixing lengths in the grid is shown in 
Fig. [2] The majority of the models have an a roughly between 
1.5 and 2.5, while we had to cut off at a — 3.5 for models with 
r eff % 3 000 K, log g % 3 and at a = 1 for those with T eS % 
6 000 K, logg £ 5. For almost all models with T eff £ 8 000 K the 
mixing length is a = 1 . 

Because all the other changes introduced with the new 
PHOENIX library, we unfortunately are unable to comment on 
the impact of our choice of this varying alpha parameter on the 
final spectra. 

2.3.3. Micro-turbulence 

Micro-turbulence is an ad-hoc parameter that is used to match 
synthetic spectra with observations. We use it in the sense that 
the large-scale (macro-) turbulent motion triggers small-scale 
(micro-) turbulent motion below the length scale of the mean 
free path of the photons. Unlike macro-tur bulent motio n, this 
motion affects the strength of a spectral line (Gray 2005), since 
it directly influences the line opacities, so it has to be included 
in the line forming process and cannot be applied afterwards. 
Following this scheme, micro-turbulence is strongly related to 
the macro-turbulent motion and we therefore use the relation 

V'micro = 0.5 • <V conv ) (7) 

to determine its amount. This empirical approximation follows 
from 3D radiative h ydrodynamic investigations of cool M stars 
dWende et al.ll2009h . but also gives reasonable values for hotter 
stars (see Fig. [3). The mean value of the macro-turbulence (v con v) 
is taken from the PHOENIX atmospheres prior to spectral syn- 
thesis. The convective velocities from all layers are averaged in a 
way that only layers with non-zero values are included and then 
divided by a factor of two, i.e. micro-turbulence is assumed to 
be half of the mean convective velocity in the photosphere. 

Since micro-turbulence mainly affects the strength and shape 
of the lines and not the structure of the atmosphere, it is only in- 
cluded in the computation of the high resolution synthetic spec- 
tra and not in the underlying model atmospheres. We did not 
attempt to further converge the atmospheres with the micro- 
turbulences derived from the previous iterations, which would 



have cost additional computational resources without noticeably 
changing the structure. Despite the lack of a self-consistent de- 
scription derived from first principles, this empirical description 
removes micro-turbulence as a free parameter. 

Figure [3] shows the distribution of micro-turbulences for 
models with solar abundances. As one can see, for cool stars 
our micro-turbulences are a lot smaller that the 2 km/s frequently 
assumed in model atmo spheres. Over-plotted are observational 
data p oints taken from Land street et al.l (120091) and Takeda Tt al.l 
(2008). The agreement with our values for stars on the upper 
main-sequence and the red giant branch is gratifying. For ef- 
fective temperatures above 6500K, however, convective tur- 
bulence predicted by our models decreases steadily and can- 
not reproduce the signi ficant broadening velocities observed by 
lLandstreet et al.1 (120091) in F and A stars. 

We note however that the lLandstreet et all (12009) study does 
not display a unique relation of photospheric line broadening 
and atmospheric parameters, and includes a number of peculiar 
stars - e.g. all main sequence stars with v m i cro > 3 km/s in then- 
sample are Am stars. They also conclude that the line shapes of 
these stars cannot be understood in a simple micro-turbulence 
concept. These findings are e ssentially confirmed by RHD sim- 
ulations (Frevtag et al.ll2012l and references therein), but these 
still sample the hottest stars with convectively unstable atmo- 
spheres only very sparsely. 

For giant stars with r e jf « 6 500 K we encountered strong 
numerical instabilities in the PHOENIX convection module as a 
consequence of to the diminishing extension of the convectively 
unstable layer and had to disable it completely. For consistency 
with the rest of our grid, we have chosen to publish these spec- 
tra with zero v^cn,; the users are advised that they cannot replace 
more detailed modelling using the micro-turbulence as a free pa- 
rameter for those stars. 

Other unexpected results appear in the region below 3 000 K, 
where we see another peak in the micro-turbulence. Those mod- 
els develop a second convective zone at the outer edge of the 
atmosphere, where the convection velocity increases towards 
lower optical depths, which causes our method to produce values 
for the micro-turbulence that are higher than expected. A simi- 
lar behaviour h as been obse r ved in 3D radiative hydrodynamic 
simulations by IWende et al.1 (120091) . Although the effect of this 
second convective zone on the spectrum is small, it does intro- 
duce changes as a result of our new self-contained method for 
calculating the micro-turbulence. Thus, we will provide an up- 
date for the spectral library for stars with T e ff ^ 3 000 K within 
a couple of months, in which we will limit the range of optical 
depths used for calculating the micro-turbulence. 

2.3.4. Equation of state 

The Astroph ysical Chemical Equilibrium Solver (ACES, 
Barman 2012) equation of state (EOS) that is used starting with 
PHOENIX version 16 is a state-of-the-art treatment of the chem- 
ical equilibrium i n a st ellar atmosphere. It uses the method of 
ISmith & Missenl £l982) with new experimental and theoretical 
thermodynamical data (Barman 2012) for 839 species (84 ele- 
ments, 289 ions, 249 molecules, 217 condensates). 

In each layer, the chemical equilibrium for all used atomic 
and molecular species is computed in dependence on the pres- 
sure, temperature, and density. After computing the atmospheric 
structure with the equations of radiation and hydrodynamics, a 
new chemical equilibrium is calculated in a second step. This 
iterative process continues until convergence. We did not not in- 
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Fig. 3. Distribution of micro-turbulence for models with solar abundances for different effective temperatures T e s and surface gravi- 
ties log g (colour-coded). We had to disable the convection for the giants with r e ff > 5 500 K and therefore no micro-turbulence was 
applied to these spectra. In and around this area we o bserve numerical instab ilities ca using single models to show differing results. 
Over-plotted are observational data points taken from Land stre et et al.l J2009) (+) and Take da et al.l ((2008) (x). 



vestigate which changes in the spectra were introduced specifi- 
cally by using the new ACES equation of state. 

2.3.5. Element abundances 

The element abundan ces are scaled solar abundances taken from 
Asplund et al. (2009), where the abundances in both the photo- 
sphere and in meteorites are given. With the exception of Li, 
which is destroyed in the solar interior, and volatile elements 
that are depleted in meteorites, both sets of abundances agree 
within the errors; apart from the cases mentioned above, we have 
therefore chosen the measurements with respective smaller un- 
certainties (see Table [H to best rep roduce the proto-solar ele- 
mental composition (Lodders 2003|). 

2.4. Vacuum wavelengths 

All our spectra are provided at vacuum wavelengths A v . dc . If air 
wavelengths A^ are required, one can use the following equation 
taken from ICiddorl (11996I) . where it is applied for A > 2 000 A 
only 



/tail- — 



with 



/ ' 



(8) 



The difference in flux between an interpolated spectrum and 
one calculated with PHOENIX is usually less than 1%, but in- 
creases for lower temperatures and can reach up to 10% for the 
coolest stars in the grid. 



3. Quality of the spectra 

Many changes have been implemented with respect to previous 
PHOENIX models. For example, the more detailed treatment of 
the chemical equilibrium in the new EOS strongly affects the 
stellar structure and results in different line and molecular band 
strengths, which can introduce significant differences in compar- 
ison to older PHOENIX model spectra, especially for M stars as 
discussed below. But major changes can also be expected from 
the new list of element abundances and the new parametrisations 
for the mixing length and the micro-turbulence. Consequentially 
we observed some significant differences between spectra from 
previous PHOENIX grids and from this one. 

We compared some temperature profiles with those of the 
original models that have been used as a starting point for the 
new library. For t | 1 they match very well and we only see 
differences for t » 1, which are irrelevant for the morphology 
of the final spectrum. 



o" 2 = 



10 4 



, 0.05792105 0.00167917 
' + 238.0185 -cr 2 + 57.362-o- 2 ' 



(9) 
(10) 



2.5. Interpolation 

A few spectra are missing in the grid where PHOENIX was un- 
able to compute the structure of the atmosphere for various rea- 
sons. Where possible, we provide interpolated spectra. 

The interpolation is done point-wise, i.e. every wavelength 
point is interpolated independently using a cubic spline created 
from neighbouring spectra. The interpolation is done along the 
r e ff axis only. 



3.1. M stars 

In Fig.|4]we show a comparison between model spectra from our 
new PHOENIX grid and from a previous one using the old EOS 
AMES version 2.6 at the region around the TiO e-band. In both 
cases, the Cond models are used. We show spectra at effective 
temperatures of T e ff - 2 900K and logg = 5.0, typical values 
for an M6 main sequence star. The TiO bands are known as ro- 
bust temperature indicators, because they only weakly depend 
on gravity but are ve ry sensitive to temperature (Mo hantv et al.l 
2004; lReinersll2005l) . In the new grid we now see TiO e-bands 
that are significantly deeper than those in older PHOENIX ver- 
sions for identical atmospheric parameters. In this example, the 
difference in temperature derived from the two models would be 
about 200 K. 
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Tabl e 4. Element abundances used in grid, taken from 
Asnl und et al.l (120091) . All the abundances in the table marked 
with an m are from meteorites, while the others are photospheric. 



z 


Elem. 


Abund. 




Z 


Elem. 


Abund 




1 


H 


12.00 




44 


Ru 


1.75 




2 


He 


10.93 




45 


Rh 


1.06 


m 


3 


Li 


3.26 


m 


46 


Pd 


1.65 


m 


4 


Be 


1.38 




47 


Ag 


1.20 


m 


5 


B 


2.79 


m 


48 


Cd 


1.71 


m 


6 


C 


8.43 




49 


In 


0.76 


m 


7 


N 


7.83 




50 


Sn 


2.04 




8 


O 


8.69 




51 


Sb 


1.01 


m 


9 


F 


4.56 




52 


Te 


2.18 


m 


10 


Ne 


7.93 




53 


I 


1.55 


m 


11 


Na 


6.24 




54 


Xc 


2.24 




12 


Mg 


7.60 




55 


Cs 


1.08 


m 


13 


Al 


6.45 




56 


Ba 


2.18 




14 


Si 


7.51 




57 


La 


1.10 




15 


P 


5.41 




58 


Cc 


1.58 




16 


S 


7.12 




59 


Pr 


0.72 




17 


CI 


5.50 




60 


Nd 


1.42 




18 


Ar 


6.40 




62 


Sm 


0.96 




19 


K 


5.08 


m 


63 


Eu 


0.52 




20 


Ca 


6.34 




64 


Gd 


1.07 




21 


Sc 


3.15 




65 


Tb 


0.30 




22 


Ti 


4.95 




66 


Dy 


1.10 




23 


V 


3.93 




67 


Ho 


0.48 




24 


Cr 


5.64 




68 


Er 


0.92 




25 


Mn 


5.43 




69 


Tm 


0.10 




26 


Fe 


7.50 




70 


Yb 


0.92 


m 


27 


Co 


4.99 




71 


Lu 


0.10 




28 


Ni 


6.22 




72 


Hf 


0.85 




29 


Cu 


4.19 




73 


Ta 


-0.12 


m 


30 


Zn 


4.56 




74 


W 


0.65 


m 


31 


Ga 


3.04 




75 


Re 


0.26 


m 


32 


Ge 


3.65 




76 


Os 


1.40 




33 


As 


2.30 


m 


77 


Ir 


1.38 




34 


Sc 


3.34 


m 


78 


Pt 


1.62 


m 


35 


Br 


2.54 


m 


79 


Au 


0.80 


m 


36 


Kr 


3.25 




80 


Hg 


1.17 


m 


37 


Rb 


2.36 


m 


81 


Tl 


0.77 


m 


38 


Sr 


2.87 




82 


Pb 


2.04 


m 


39 


Y 


2.21 




83 


Bi 


0.65 


m 


40 


Zr 


2.58 




90 


Th 


0.06 


m 


41 


Nb 


1.46 




92 


U 


-0.54 


m 


42 


Mo 


1.88 













3.2. Comparison with Kurucz models 

In Figs. [5] and [6] the PHOENIX synthetic spectra are com- 
pared to corresponding ATLAS9 spectra (ICasteili & Kurucz] 
120031) f or two different effective temperatures. Note that they 
used the lGrevesse & Sauvall(ll998l) solar abundances, wh ile we 
used the more recent abundances bv lAsplund et al.l d2009). This 
will introduce systematic differences between the spectra. The 
PHOENIX spectra have been convolved with a Gaussian ker- 
nel with a FWHM of 20 A for wavelengths smaller than 1 /an 
and with a Gaussian with a FWHM of 50 A for wavelengths 
larger than 1 //m in order to match the spectral resolution of the 
ATLAS9 spectra better. In general the compared models agree 
well. For the 8000 K spectra we observe a ~3% lower flux at 
every single wavelength point in the PHOENIX spectrum than 
in the ATLAS9 spectrum, which is therefore just a difference in 
total luminosity that can easily be explained by the varying mass 
in our new library. 



ACES-AGSS-COND 

AMES-cond.v2.6 




0.4 I i i I I I i i 

8420 8440 8460 8480 8500 8520 8540 

Wavelength [A] 



Fig. 4. New PHOENIX ACES model spectra (full line) com- 
pared to PHOENIX AMES v2.6 spectra (dotted line); both mod- 
els are calculated for a typical M6 main sequence star (T e ff = 
2 900K, logg = 5.0) and smoothed to an effective resolution 
of R = 10000. The ACES model spectra produce significantly 
deeper TiO bands than earlier PHOENIX models. 




t 1 1 1 1 1 1 

5000 10000 15000 20000 25000 30000 

Wavelength [A] 

Fig. 5. Comparison between PHOENIX (full line) and ATLAS9 
(dotted line) model spectra for T e s = 4000 K, log g = 4.0 and 
solar abundances. The PHOENIX spectrum has been convolved 
with a Gaussian kernel with a FWHM of 20 A for wavelengths 
smaller than 1 /mi and with a Gaussian with a FWHM of 50 A 
for wavelengths larger than 1 /zm. 



Table 5. Input parameters for PHOENIX and ATLAS9 models 
for Sun and Vega. Please note tha t \Fe/H] in the PHO ENIX 
models is give n by abundances from Asplun d et al.ld2009l) . while 
ATLAS9 uses lAnders & Grevessel(ll989l) . 



Parameter 


Sun 




Vef 


;a 




PHOENIX 


ATLAS9 


PHOENIX 


ATLAS 9 


T t s 


5778 K 


5777 K 


9550 K 


9550 K 


logg 


4.44 


4.44 


3.95 


3.95 


[Fe/H] 


0.0 


0.0 


-0.3 


-0.3 
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Comparison of spectra for T eff = 8000K and Ic 



Table 6. Comparison of colour indices for different solar spectra. 

















if 


































1 






























i 


PHOENIX 

ATLAS9 
l 









5000 5500 
Wavelength [A] 



Fig. 6. Comparison between PHOENIX (full line) and ATLAS9 
(dotted line) model spectra for = 8000 K, logg = 4.0 and 
solar abundances. The PHOENIX spectrum has been convolved 
with a Gaussian kernel with a FWHM of 20 A. 



H a line in solar spectrum 




6545 6550 6555 6560 6565 6570 6575 6580 6585 
Wavelength [A] 

Fig. 7. Comparison between PHOENIX (full line) and ATLAS9 
(dotted line) model spectra for the Sun in full resolution (R = 
500 000) in the upper plot. Below the PHOENIX spectrum (full 
line) is compared to a solar spectrum observed bv lWallace et alJ 
(2011. 



3.3. Sun 

In the upper plot of Fig. [7] a comparison between a PHOENIX 
and an ATLASS0 solar spectrum in full resolution is shown for 
the region around H a . Table[5]lists the input parameters for both. 
The overall agreement is good, but the line profile is slightly 
different: in the PHOENIX spectrum the line is not as deep as in 
the ATLAS9 spectrum, which is compensated by broader wings. 
The lower plot of Fig. [7] com pares the P HOENIX spectrum to a 
solar spectrum observed by Walla ce et al.l (1201 lb . Obviously the 
real line is even broader and thus the PHOENIX spectrum is a 
better match than the ATLAS9 one. 



Spectrum 


U-B 


B-V 


V-R 


R-I 


V-I 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


ATLAS9 


0.069 


0.654 


0.443 


0.324 


0.767 


PHOENIX 


0.065 


0.633 


0.441 


0.327 


0.768 


CALSPEC 


0.083 


0.645 


0.446 


0.325 


0.771 




6000 7000 
Wavelength [A] 



Fig. 8. Comparison between PHOENIX (full line) and ATLAS9 
(dotted line) model spectra for Vega, both convolved to 20 A 
FWHM. 



Unfortunately, these comparisons also show one of the major 
drawbacks of synthetic spectra, i.e. they are missing some of the 
absorption lines present in the observed spectrum, e.g. the line 
at 657 1 A, which is present in the ATLAS9 spectrum but not in 
the PHOENIX spectrum. 

Some colour indices calculated from the ATLAS9 and the 
PHOENIX spectrum are listed in Table [6] In addition, the 
colours for the solar spectrum from CALSPEC, which is a 
set o f composite stellar spectra used for HST c alibrations (see 
e.g. lTumshek et alJll990tlColina & Bohlrnll 19971) are given. The 
transmission cu rves for the used filters have been adapted from 
Uohnsonl(ll965lPl 

Altogether our solar spectrum seems to be quite accurate, 
which is gratifying for an all-purpose spectral library without 
special optimisations to match the solar spectrum well. 



3.4. Vega 

In Fig. [8] we show a comparison between spectra for Vega, cal- 
culated both with PHOENIX and ATLAS'S and the input pa- 
rameters given in Table [5] The PHOENIX model has been in- 
terpolated from the existing grid instead of calculating a new 
one. Since the mass of that model differs from the real Vega be- 
cause of our parametrisation of the stellar mass, the flux was 
scaled to match that of the ATLAS9 model. Further discrep- 
ancies are introduced by the differences in helium abundance, 
micro-turbulence, etc. Nevertheless, the overall shapes of the 
spectra match well. 

Table|7]lists some colour indices. Since the magnitudes have 
been calculated in the VEGAmag system with the ATLAS9 
model as reference, all colours are zero for these magnitudes. 



http : //kurucz . harvard . edu/sun . html 



3 http : //obswww. unige. ch/gcpd/filters/fil08 .html 

4 http : //kurucz . harvard. edu/stars/vega/ 
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Table 7. Comparison of colour indices for different spectra of 
Vega. 



Spectrum 


U-B 


B-V 


V-R 


R-I 


V-I 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


ATLAS 9 


0.000 


0.000 


0.000 


0.000 


0.000 


PHOENIX 


-0.019 


0.005 


0.007 


-0.002 


0.005 


CALSPEC 


0.027 


0.006 


0.006 


0.008 


0.014 



The results from our interpolated PHOENIX spectrum are in 
good agreement. 

4. FITS files 

4.1. Download 

All the spectra presented in this paper are available for down- 
load from our homepage^ Single high-resolution spectra can be 
downloaded using a web form by specifying temperature, sur- 
face gravity, metallicity and alpha element abundance. In order 
to keep the amount of storage space to a minimum, these files 
contain only the flux; the wavelength grid must be obtained sep- 
arately once for all spectra. 

Furthermore we provide archives with full sub-grids of the 
library convolved to more convenient resolutions: 

1 . AA = 1 A in the optical wavelength range from A = 3 000 A 
to 10 000 A. 

2. R - 10000, i .e. the resolution of the X-Shooter Spectral 
Library XSL (IChen et all 1201 ll) . with the full wavelength 
range of X-Shooter, namely A = 3 000 A to 24 800 A. 

These spectra are over-sampled by a factor of ten, so the 1 A grid 
has a sampling rate of 0. 1 A, for example. 

4.2. Naming scheme 

The spectra are provided in a FITS file with the following nam- 
ing scheme: 

<grid>/<subgrid>/ 
(n) lte<Tef f xlog (g) ><subgrid> . <grid>-HiRes . fits 

where <grid> is the name of the grid. As for the 
spectra described in this paper, this will always be 
PHOENIX-ACES-AGSS-COND-2011. 

The <subgrid> describes all the parameters in addition to 
r e ff and logg and is empty when there are none. An example 
would be . Alpha=+0 . 30 for an alpha element enhancement of 
+0.3 dex. Abundances equal to those of the Sun are never in- 
cluded, i.e. there is no . Alpha=+0 . 00. Be aware that [Fe/H] is 
given as first parameter in the name of the sub-grid, not Z. 

The filename itself always starts with either lte or nlte, 
describing whether the atmosphere has been calculated in LTE 
or NLTE. 

The effective temperature <Te££> is always given by five 
digits with leading zeros, if necessary. The surface gravity 
<log(g)> is denoted by its negative with two digits after the 
decimal point. In the filename itself the Z describing the metal- 
licity is dropped, since the metallicity is explicitly given for all 
files. 

Two examples for the naming scheme would be 
5 http : //phoenix . astro . physik . uni - goettingen . de/ 



* PHOENIX- ACES- AGSS-COND- 20 1 1/Z- 1 . . Alpha=+0 . 30/ 

lte08000-2 .00-1.0. Alpha=+0 . 30 . 
PHOENIX-ACES-AGSS-COND-2011-HiRes.£its 

* PHOENIX-ACES-AGSS-COND-2011/Z-0 . 0/ 

lte©68©0-4. 50-0.0. 

PHOENIX-ACES-AGSS-COND-2011-HiRes.£its 

4.3. File content 

The files always contain one single primary extension, which 
holds the flux of the spectrum in units of [erg /s /cm 2 /cm] on the 
stellar surface. The files for the high-resolution spectra contain 
a FITS keyword WAVE, which holds a reference to the file that 
stores the wavelength grid and is saved in the same format as the 
spectra, but with wavelength points instead of a flux array. 

For the medium-resolution spectra, the wavelength grid is 
provided by a set of four FITS keywords. CRVAL1 defines the 
wavelength at pixel CRPIX1, which is always 1 for our spectra. 
The step size is given by CDELT1 and the number of points 
by NAXIS1. The type of the wavelength grid is defined by 
the CTYPE1 keyword — for this grid it will always be one of 
the two possible values for vacuum wavelengths, i.e. WAVE or 
WAVE-LOG, where the latter denotes a logarithmic scale. With 
this information the whole wavelength grid can be calculated 
easily. 

4.4. Additional FITS keywords 

In addition to the default FITS keywords, we added some more 
for documenting the atmospheric parameters for the spectrum. 
The basic parameters are given as 

- PHXTEFF: [K] effective temperature, 

- PHXLOGG: [cm/s 2 ] log (surface gravity), 

- PHXM_H: [M/H] metallicity (rel. sol. - Asplund et al 2009), 

- PHXALPHA: [a/M] alpha element enhancement. 

Additional parameters for the atmosphere are 

- PHXDUST: [T/F] Dust in atmosphere, always F for this 
grid, 

- PHXXI_L: [km/s] micro-turbulent velocity for LTE lines, 

- PHXXI_M: [km/s] micro-turbulent velocity for molec lines, 

- PHXXI_N: [km/s] micro-turbulent velocity for NLTE lines, 

- PHXMASS: [kg] Stellar mass, 

- PHXREFF: [cm] Effective stellar radius, 

- PHXLUM: [W] Stellar luminosity, 

- PHXMXLEN: Mixing length. 

Furthermore, we included some information about the 
PHOENIX run: 

- PHXBUILD: Build date of Phoenix, 

- PHXVER: Phoenix version, 

- DATE: [local] finishing date of model calculation, 

- PHXEOS: Equation of state. 

As mentioned before, there are models within the grid for which 
we manually disabled the convection. For those we added an- 
other keyword PHXCONV with the value F. For the high- 
resolution spectra there is also a second FITS extension in the 
file containing a binary table listing all the element abundances 
that have been used in the model. 

The files containing interpolated spectra only include the 
date of their creation and those FITS keywords from the lists 
above that specify the basic stellar parameters (T e ff, logg, 
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[Fe/H], and [a/Fe]). Furthermore, for these spectra we added 
the keyword INTERPOL with the value T. 



5. Summary 

The presented PHOENIX grid contains a comprehensive set 
of synthetic spectra that allows a detailed analysis of observed 
spectra for a wide range of applications. An extension taking into 
account NLTE treatment of important elements up to 25 000 K is 
a work in progress. Both resolution and wavelength range match 
existing and upcoming state-of-the-art instruments. 
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